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We investigate the effects of WW rescattering through strong anomalous four-vector boson cou- 
plings. In the I — 1, J — 1 channel, we find a resonance with a mass of approximately 200GeV 
and a width of less than 12GeV. In an application to pion physics we find a small correction to the 
KSRF relation. 

I. INTRODUCTION 

Within the standard model of electroweak interactions the gauge principle fully determines the selfcouplings of 
the vector bosons. Therefore the measurement of these couplings is of prime importance, as deviations would be 
an indication of new physics. The gauge principle predicts well-defined three and four-vector boson selfcouplings. 
Deviations of these values can be described in a gauge-invariant way in the Stiickelberg formalism Within the 
Stiickelberg formalism the standard model is described as a gauged nonlinear sigma model. This implicitly assumes 
that the Higgs particle does not play a fundamental role. Alternatively one can see it as the m# — * oo limit of the 
standard model. If indeed anomalous vector boson couplings are present, this should be a reasonable assumption, 
since in that case strong interactions should be present. The triple vector boson couplings are severely constrained 
by the measurements at the CERN e + e~ collider LEP-200 and the Fermilab Tevatron 0]- Also, indirect limits from 
LEP-100 and (g — 2) M exist. Altogether experiments indicate they should be small. This is not too surprising, as it is 
very hard to construct a model that would give rise to large effects. The reason is that within the three-vector boson 
couplings there is always an interplay between longitudinal and transversal vector bosons. 

For the four- vector boson couplings the situation is somewhat different. Here one can write down vertices that 
contain longitudinal vector bosons only. These vertices correspond to the Goldstone boson sector of the theory and it 
is much easier to generate strong interactions in this sector. Such results come typically through intermediate heavy 
Higgs boson exchange. An example of such a model is given in ||, where the strong interactions are generated via 
singlet effects in the Higgs sector. Also, in the standard model the two-loop heavy Higgs correction in the four-vector 
boson j| couplings is an order of magnitude larger than in the three- vector boson couplings || . About the four- vector 
boson couplings much less is known than about the three- vector boson couplings. Direct experiments probing these 
interactions do not exist at present. They can at the moment only be tested through radiative corrections in the p 
parameter. These corrections can be calculated within perturbation theory with a cutoff A . Within the four-vector 
boson couplings one should distinguish between two types. In the standard model there is an extra global SU_r(2) 
symmetry when the hypercharge is turned off. For the anomalous couplings this is not necessarily the case. The 
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couplings that violate SUfl(2) even in the absence of hypercharge give quartically divergent corrections to Sp 
and should therefore be negligibly small. Physically, this means that the underlying strong interactions, generating 
the anomalous couplings, should respect the SUr(2) invariance. This leaves only two operators that preserve SU#(2) 
invariance in the absence of hypercharge. If one could use simple cutoff perturbation theory these operators still 
give quadratically divergent corrections to Sp |(|f7j and are sufficiently suppressed to give only small effects in future 
colliders. In [0 the quadratic divergences were ignored and therefore the limits are weak. 

This leaves only the possibility that the anomalous couplings are so large that perturbation theory cannot be trusted 
and therefore the low energy limits are invalidated. It is precisely this case that we study in this paper. We assume 
that anomalous couplings are present which preserve the SUfl(2) symmetry of the theory and assume the coefficients 
for the corresponding operators to be large. To simplify the calculation we actually ignore hypercharge altogether and 
work in the SUz,(2)xSU#(2) model. This leads to a simplification because we can study different isospin channels 
separately. Because the interactions are assumed to be strong, vector boson scattering cannot be described by the 
tree-level vertex. We therefore perform a resummation of loop graphs. In most channels we find no particularly 
interesting effect. The only exception is the I = 1, J = 1 channel, where a resonance is found. In the case of very 
strong anomalous couplings, which we assume, the resonance can be quite close to the two- vector boson threshold. The 
coupling of the resonance to the vector bosons is found to be suppressed by the cutoff and could be small. Dependent 
on the parameters the resonance could be visible even at LEP-200. 

The paper is organized as follows. In section || we present the model. In section [II we perform the calculation of 
the bubble sum. In section IV we discuss the results. The appendices contain technical details. 



II. THE MODEL 



A. Lagrangian 

We work in a pure massive SU(2) gauge theory and introduce the anomalous couplings in a gauge-invariant way 



using the Stuckelberg formalism [jy. That is, we write the theory as a gauged nonlinear sigma model: 

C„ = -±Tr(WV^n + ml Tv(V^Vn (1) 

with 

W, ee \r a W« = W\ , (2) 

W MV ee \r a W; v = d„W v - DvW^ + ig[W„ W v ] = Wl v , (3) 

V,ee-1(D,U)W = Vl, (4) 

D^UEEd^U + igW^U, (5) 

U = exp(iv a T a ) (6) 

with real fields v a , i.e., the SU(2)-valued field U describes the Goldstone degrees of freedom. m WQ is the mass of the 
three gauge bosons in absence of the higher covariant derivative terms to be added below. The anomalous couplings 
are then introduced as 

Cano = g4(Tr[V„V v }) 2 + <7 5 (Tr[VV 

= g 4 (TrlV„Vt}) 2 + g 5 (Tr[V»Vl ]) 2 

= ff 45- 4 {Tr[(^f/)(L>,f/)t]} 2 + g 5 g-HTr{(D»U)(D,Uy}y . (7) 



2 



Their form is determined by the requirement that they conserve CP, are not accompanied by three-vector boson 
couplings (on which the limits are much more stringent, as was mentioned in the introduction) and that they are 
invariant under the custodial SUr(2) symmetry U — ► UUr with U £ SU(2). 

To regulate higher-than-logarithmic divergences, we introduce higher covariant derivative terms through 

1 777 2 

C hcd = Tr[(D a W^(D a Wn] - ^Tr[(D a V,)(D a Vn] (8) 

with 

D a W^ = d a W^ + ig [W a , W^} , (9) 

Aw and Ay effectively being momentum-cutoffs. These are the unique dimension 6 higher derivative propagator terms 
and are further discussed in 0. 

We work in unitary gauge, where we have U = 1, v a = 0, and therefore 



The signature of our metric is 



v M = wv (io) 

w = diag(l, -1,-1,-1). (11) 



B. Feynman Rules 

The W propagator in unitary gauge with higher covariant derivatives is 

A^(k) = A? r (k 2 ) (g, v + A,T(fc 2 )%^ (12) 



with 



K(k 2 ) - (fc2 _ m2 _ )( ^ 2 _ o j , (13) 



m wo k 



A 2 



(14) 




A 2 



The Feynman rule for the anomalous four-vector boson couplings is 

■c,7 



O(A ) 

4o + o(a- 2 ) 



A 2 



(15) 



a, a 



b,/3 




= i{ S ab S cd [2g 5 g a! 3g^s + gA(g ai g/35 + 9a59f3~f)] 

d, 6 +SacSbd[^959a-y9l3S + .94 (.9 

+5 ad S bc [2g 5 g aS g /3l + 9i{g ai gi3& + 9af39 7 s)]} 

Trabcd 

= V af3y6 • 



(16) 
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III. BUBBLE SUM 
A. Definitions 



Assuming the anomalous couplings to dominate the gauge coupling, we compute the "bubble sum" 



a, a 



b,(3 




c, 7 a, a 

d, 5 b,(3 

a, a 

b, (3 





c, "f a, a 

+ 

d, 6 b,/3 

c, 7 a, a 
+ 

d, 6 b,/3 





c, 7 a, a 

+ 

d, 5 b,P 
c, 7 a, a 



d,6 b 7 (3 





c, 7 

d, S 

c, 7 

d, 5 



d(0) abed j-y(l) abed jj(2) abed 
seabed 



(17) 



where all lines represent gauge bosons and where the vertices contain only the anomalous part. Note that (|l7j ) is 
invariant under each of 



(a, a) <-> (&,/3) , 
(0,7) <-> (d,<5) , 
(a,a,b,/3) <-> (c,7,d,<5) . 



(18) 
(19) 
(20) 



To express equation (|]_7j) in a purely algebraic form, we split each additional bubble into a vertex part ( |16| ) and a 
propagator part 



b,0 



c,7 



d,<5 



a, a 



b,(3 



c, 7 a, a 



d, 5 b,(3 




c,7 



d,5y 



(21) 



where we have imposed the relevant symmetries (|18|)-(|2C|). equation (|l7j) can now be written as 

a,a . . c,7 

&,j9 / \ d,S 

a, a 



b,P 




c, 7 a, a 



d,<5 




A 




c, 7 a, a 



d,<5 




A 





c,7 



d,<5 



(22) 



Since the vertex is independent of any momenta, the integration over the loop momenta can be done in the propagator 
part. Abbreviate 



d d p 
{2n) d 



(23) 



with the renormalization scale /i and use dimensional rcgularization together with the modified minimal subtraction 
(MS) scheme throughout. Define 
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p 



abed _ 1 



b,0 



bil 



A 



. c,7 

V)(. 



</a 7 p 2 



A 



it" 



((P + fc) 2 )^ 



Ig \f l p 
(p+k) f i(p+k) s \ , A 



2^1 (p+fc)/a(p+fc)i 
5 



+8adSbc 



+ A%(p 2 )^ 



be 



Agr((p + fc) 2 ) ( 5 , 7 - (p+ g: ( g fc) - ) + a^((p + fc) 2 ) (p+ g: ( fc y K ] } 

Ai(fc 2 )g aT Sf 1 a5 + A 2 {k 2 ){g a pg 1 s + g a sgp~<) + A z (k 2 )(g ai kpkg + gpsk a k 7 )/k 2 
+A 4 (fc 2 )(g a/ 3fc 7 fc5 + g^skakp + g a skpk 7 + gffryk a ks)/k 2 + A 5 (k 2 )k a kffk 1 ks/k 4: 
A\{k )g a sgp 1 + A 2 {k 2 ){g a pg 1 s + g a7 gps) + A 3 (fe )(g a skpk 7 + gp~fk a ks)/k 2 
+A i (k 2 )(g a pk 7 k s + g 1 sk a kp + g aj kpk s + gpsk a k y )/k 2 + A 5 (k 2 )k a k f} k~ / k s /k 4 j , 



(24) 



where an additional factor of 1/2 has been introduced to account for the implicit symmetry factors in the bubble sum 
(|l7l), which can now be written as 



\^abcd \7-abcd . t/ a b rn n tj m n m. ' n ' rr rn ' n' cd . Trabmnprnnrn'n / jrrn / n / rn / 'n' / pm" n" m!" n'" irm"' n //; cd ■ 



(25) 



The strategy for computing Ai through A 5 is given in the appendix. Keeping only at least quadratically divergent, 
i.e. 0(A 2 ) terms for the total Ai through A5, and O(A ) terms for their real parts, the results are 



A 1 (k 2 ) = 



iA v 



96(4tt) 2 ?74 



(i 



A 2 
ln Aw 



7iA 2 v 



k 2 



ReAi(fc^) 



48(4 7 r) 2 (A 2 v -A 2 / )m 2 V0 A 2 , 1152(47r) 2 m 2 , \m 2 H 
( k 2 



+ O(A ) , 



2 ^ _ TTy/1 - 4m 2 wo /k 2 



60(4tt) 



\4m 2 v 



/ k 2 
\Am 2 v 



+ 0{K- 2 ). 



(26) 



(27) 



A 2 (fc 2 



iA v 



96(47r) 2 m 4 



wo l £ 6 
zAf^A 2 , 



A 2 

— ln^ 



7*A 2 , 



48(4n) 2 (A 2 w - A 2 v )m 2 W0 A 2 , 1152(47r) 2 m 2 



wo \""wo 



+ O(A ), 



ReA 2 (fc 2 ) 



iry/1 - Am 2 lv Jk 2 



60(4tt) : 



A 3 (fc 2 ) 



IHAy 



V4m 2 , 
k 2 



1 2 



0(A- 



576(47r) 2 m 2 , n \ m 2 



wo \ ""wo 



■O(A ) 



ReA 3 (fc 2 ) = 



ny/1 - 4m^ /fc 2 
60(4tt) 2 



/ k 2 



/ fc 2 

-13 - 



\ 4m w 



+ 0(A- 2 ), 



(28) 
(29) 
(30) 
(31) 
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A 4 (fc 2 



7iA 2 y 



576(47r) 2 m 2 V0 \m 2 V0 



O(A ) 



(32) 



ReA 4 (r) = 



2 _ Ky/l - 4m 2 , Q /fc 2 



60(4tt) 2 



k 2 



Ho 



- 1 



A 5 (fc 2 ) = 0(A°), 



+ G(A- 2 ), 



(33) 
(34) 



ReA 5 (fc 2 



tt^I - Am 2 w Jk 2 
60(4tt) 2 



/ fc 2 



/ k 2 



+ C(A- 2 ). 



(35) 



where e is defined by d = 4 — 2e with spacetime dimension d. Note that Ai(fc 2 ) and A2(/c 2 ) have the same quartically 
divergent part 



A (4) = _ 



96(47r) 2 m^ 



1 - + 5 --m^ 



e 6 



(36) 



Note furth er tha t A^(k 2 ), A 4 (fc 2 ) are only quadratically divergent and that As(fc 2 ) is at most logarithmically divergent. 
In section IIIC we will need the following linear combinations which are only quadratically divergent: 



and 



M-A,(f)l = - i , ( ^f^ +0(A ° ) 



-^(fc 2 ) - A 2 (fc 2 ) + A 3 (fc 2 ) - A 4 (fc 2 )] 



(37) 



A 2 



8(4^) 2 mt 
+0(0,-2), 

where the shorthand notation 



k 2_ 4m 2 In(A^/A^) , z^(P + 4m 2 J(fc 2 -4m 2 , ) 3 / 2 



W °1-A 2 V /A 2 W 



3VtfA 2 v 



0(m, n) = ReO(A m ) + ImO(A") 



(38) 
(39) 



has been used. 



B. Tensor Structure of and P$ff s 

Now let us analyze the tensor structures that can appear in the bubble sum. The SU(2) tensors in ( |l6|) and ( |24| ) 
are 

SabScd , S ac Sbd , S a d$bc , (40) 

which can be rewritten into the following linear combinations corresponding to isospin 0, 1 and 2 contributions 
so, s\ and S2, respectively, in the (a, b) <-> (c, d) channel, where we also indicate the parity under the symmetry 
transformations (p^)-([2C|). 





a <-> b c <-> d (a, b) (c, d) 


«0 = 

51 = 

52 = 


jfiabScd 

2^acSbd — Sadfoc) 
— \$ab$cd + \ {5ac$bd 


+ SadSbc) 


+ + + 
+ 

+ + + 



(41) 
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These relations can be inverted to give 

SabScd = 3s , (42) 

SacSbd = s + Si + s 2 , (43) 

SadSbc = s - Si + s 2 . (44) 

Define the product s^si by (sk) a bmn(si) mnc d- This can be analyzed in terms of the Sk again such that we have an 
algebra 

SfcS; = OklrnSm . (45) 

It is easy to see that the akim are given by 

f 1 k = l = m, 

aklm = { otherwise, (46) 

which merely means that the different isospin channels do not mix. This can be illustrated by defining the matrix S 
with elements Ski = SkSi, 

' s 

S=\ si | . .47) 

S2 

where empty entries are vanishing, and observing the absence of non-zero off-diagonal elements. 
Going to an isospin basis 

\±)=W^ = ^(W^TiW^), (48) 
|0>=Z M = iW*, (49) 

we can use Clebsch-Gordan coefficients to write 

So = |0,0)(0,0|, (50) 
i 

Sl =J2 IMXM, (si) 

fc=-l 

2 

,s 2 =]T \2,k){2,k\, (52) 

fc=-2 

where the first entry means the total isospin and the second its three-component. 
The Lorentz tensors in (|l6|) and (|24|) are 

g a f3g~/8 , g ai 9fi& > ga&g^ > (53) 

gapkyks , gajkpks , g a skpk^ , g^sk a kp , gpskaky , gpjkaks , (54) 

Define 

C = <V " k ^/ k " ■ (56) 
Then the tensors in (|5^)-((5^) can be expressed in terms of the linear combinations 
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1 n tr n tr 
d-l9a(39 7 S ' 



d-lk 2 ' 
£4 = 7^7 k a kp k*y ks 1 

*5 = --^19^139^8 + k^a-ygjfs + 9aS9%) 1 

j- 1 /^.ir „tr „tr „tr \ 

f 6 — 2(9aj9pS - 9aS9f3y) ; 

£7 = 2p'(5a'7^/3^(5 + g^S^akj + g^gk^kj + S/^^a^) j 

is = ^i{g%k p k s - gfighak-f - g^kpk^ + gf 7 k a k s ) , 
2w{9aykf3ks — g t ^k a k 1 + g t ^k^k 1 — g^kakg) , 
zkzigorykpks + g^gkakj — g^gkfjk^ — g^kaks) ■ 



t\0 



Define the product tkU by (tk) a p^ (iij^^j and get the algebra 
The only non-zero Tkim are 

Tl,l,l = r l,2,2 = 72,3,1 = r 2,4,2 = 73,1,3 = T 3,2,4 = 74,3,3 = 74,4,4 
= 75,5,5 = 76,6,6 

= 77,7,7 = 7-7,8,8 = 7-6,9,7 = 7-8,10,8 = 7-9,7,9 = 7-9,8,10 = 7"10,9,9 = Tio,lo,10 = 1. 

To make the structure of this algebra more transparent, define the matrix T with elements T^i = tkU, 

\ 



T = 



(57) 
(58) 
(59) 
(60) 
(61) 
(62) 
(63) 
(64) 
(65) 
(66) 

(67) 



(68) 




(69) 



where all empty entries are vanishing and where the lines have been drawn to guide the eye. 

Since the different isospin channels decouple from each other, let us decompose each relevant tensor 1°^ 



e.g. 



■^a^yS = X = X S + XiSi + X 2 S2 ■ 



V = V s + V 1 s 1 + V 2 S2, 
P = Paso + Pi si + P 2 s 2 , 
B^=B^ Sa + B^ Sl + B^ S2 . 



by writing 
(70) 



(71) 
(72) 

(73) 



The Xi can be decomposed into the tk by writing 



Xi — Xj)tk , 



(74) 



where k runs from 1 through 10. For example, we have for the anomalous couplings 
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Vq = 2i{[(d+3)g 4 + (3d-l) fl5 ]ti + Vd=l{g A + 3g 5 ){t 2 + ta) + 5( 5 4 + Qz)U + 2(2g 4 + fl B )(*B + *?)]} ■ 



Vi =-2i (ff4-2fl 6 )(*6+tio), 

V^ 2 =i [2(rf 34 + 2ff B )*i + 2Vd-l gi (t 2 + t 3 ) + 4(.g 4 + g 5 )U + 2(<? 4 + 2g 5 ){t 

and for the integrated propagator part 

P = P 2 = 2(Ax + dA 2 )t 1 + 2Vd~" T(A 2 
+2(Ai + A 2 )t 5 + 2(Ai+A 2 - 
Pi = 2(Ai - A 2 )t 6 + 2(Ai - A 2 

Now we can write 

>(£+!) _ 



A 4 )(t 2 + i 3 ) + 2(A 1 
+ A 3 + A 4 )i 7 , 
A 3 - A 4 )ti • 



tr)] 



2A 2 + 2A 3 + 4A 4 + A 5 )i 4 



B\ 



B^PiVi + ViPiB^ 



(75) 
(76) 
(77) 



(78) 
(79) 



(80) 



Note that the symmetry (20) prevents fg and tg from appearing in the Vi and Pi and that therefore fl69| ) tells us that 
they will not be generated at any point in our calculation. Equation (120) also restricts t 2 and i 3 to appear only in the 



combination £ 2 + i 3 in Vi, Pi and b\ L> , which is not the case for the single terms on the right-hand side of ( |S0| ) 
Let us introduce a basis of Lorentz tensors which includes only those necessary to describe the Vj , Pj and Bj 

a <-> (3 7 «-> <5 



(i) 



1 n tr 

/ rf _ lfc a (9apk 7 ks + gjskakp) 
■jprk a kpk 7 ks + + + (81) 

— d=i9ap9ys + 2(9^9^35 + 9aS9%) 
2(9^9 f)& — 9aS9/3-y) 

(9 a~f kp ks +5/35 k a fc 7 +g^s kp k 7 +gp 1 k a ks ) 
W r (9a-y^^s+9pskock 1 —gJ s kpk 1 —g^ y k a ks) 

As in (f4l|), we have indicated the parity under the symmetry transformations (|l8|)-(^0|). 

To give the tensors u\ through u 7 a physical interpretation, let us consider them in the rest system of k^, i.e. where 



u% = 


tl 


u 2 = 




u 3 = 


u 


U4, = 




U 5 = 


te 


Uq = 


tr 


U 7 = 


tio 



(1,0,0,0) 



(82) 



and 





10 

10 

1 



(83) 



Then the Lorentz index refers to a spin-0 particle and the other three components to a spin-1 particle. In a general 
Lorentz frame, contains a spin-1 field with k^W^ = and a spin-0 field k^W^. Write the tensors in a bra and ket 
notation such that 



k a =\s)i 
kf3=\s) 2 
k 



ks 



2(S\. 



(84) 
(85) 
(86) 
(87) 



(the indices enumerate the particles) where "s" refers to a "scalar," and |1), |2), |3) are states with definite spin-up 
components in the x, y and z direction, respectively. Going to a spin basis with \s) and 
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|±l) = ^[|l) T i|2)], (88) 
|0)=*|3), (89) 



we can use Clebsch-Gordan coefficients to write 



ui=|0,0><0,0|, (90) 
U2=-(|0,0)i<*| 2 <«| + |*)i|a)2<0,0|) ) (91) 

U3 = |s)l|s>2 1<S|2<S| , (92) 
2 

u 4 =J2 \2,k)(2,k\, (93) 

fe=-2 
1 

u 5 =J2 |l,*:><l,fc| , (94) 
fe=-i 

i 

U S = ~\ E G*>ll«>2 + l*>l|*>2)(l<fc|2<*| + l(s| 2 (fc|) , (95) 
fe=-l 
1 

«7=-| E d*>il s >2 - k)i|fc) 2 )(i(fc| 2 ( S | - i(s| 2 (fc|) . (96) 
fe=-l 

where the first entry means the total isospin and the second its three-component, when two entries are present. This 
allows for the interpretation of u\, u$ and u 4 as channels for spin-0, 1 and 2 combinations, respectively, from two 
spin-1 particles. Then w 3 is interpreted as a spin-0 combination of two scalar particles and u 2 as a mixing channel 
between the spin-0 combination of two spin-1 particles and the spin-0 combination of two scalar particles, uq and u-j 
are different spin-1 combinations of a scalar and a spin-1 particle. 
Now we can write, for the vertex, 

V = 2i{[{d + 3)g 4 + (3d - l)s 5 ]«i + Vd- 1(34 + ^9b)u 2 + 5(.g 4 + .95)^3 + 2(2. 94 + .95) («4 + «e)] , 

(97) 

V x = -2i (54 - 2 35 ) («b + «7) , (98) 
V 2 = t [2(d 94 + 2g 5 ) Ul + 2^d~\g A u 2 + 4(.g 4 + g 5 )u 3 + 2(g 4 + 2g 5 )(u 4 + u 6 )] (99) 

and, for the integrated propagator part, 



P = P 2 = 2(Ai + dA 2 )u 1 + 2 v / d~^I(A 2 + A 4 )u 2 + 2(Ai + 2A 2 + 2A 3 + 4A 4 + A 5 )u 3 

+2(Ai + A 2 )w 4 + 2(Ai + A 2 + A 3 + A 4 )u 6 , (100) 
Pi = 2(Ai - A 2 )m 5 + 2(A a - A 2 + A 3 + A 4 )u 7 . (101) 



C. Transfer Matrices and Their Eigenvalues 

Writing 

B? ] = B%>u k , (102) 

we can define transfer matrices Mj, such that 

Bf k +1) = M m B<tf . (103) 

Writing 

Vi = Vi k u k (104) 
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and using that 



we can write 



B 



(0) 



Vi 



B 



Ik 



{Mf) kl V n . 



(105) 
(106) 



This can be simplified further. From (j69|) and Q81fl and the absence of tg and tg follows that m, uk, uq and wj propagate 
independently. The symmetries (18) and ([l9]) together with the parity properties noted in ([y]) and ( |si| ) necessitate 
that u§ and M7 appear only in the isospin-1 channel, while the other u k appear only in the isospin-0 and isospin-2 
channels, which we summarize in the following table. 








isospin 
1 


2 





S Ml, S U 2 , S U 3 




S2U1, s 2 u 2 , s 2 u 3 


spin 1 


s u e 


S1U5, S1M7 


S2U6 


2 


SO "4 




S 2 U4 



Even though S1M5 and S1U7 carry the same spin and isospin assignments, they do not mix. We can write 

3 



and then 



V 



fc=i 



V 1 =V 1 {5) u 5 + V^ 7) u 7 . 

3 

,(123) 



V 2 =J2V 2 r ) u k + V^u i + V 2 ^u 6 , 



k=l 



(107) 



(108) 

(109) 
(110) 



with 



B, 



(L) 
0k 



(123)Lx 



(L) 
Ik 

(L) 
24 

(L) 
26 



)klV ) 



(123) 





1=1 


R (£) 

-#04 


\(4)i T/ (4) 


-°06 


= A< 6) L vf } 


R (i) 


= A^V^ 




= A< 7) L VI 7) 
3 



(123)Lx T/ (123) 
2 Mi 



Z=l 

A 2 4) V 2 (4) 



>(6) _ 



fc = 1,2,3, 



1,2,3, 



8i (2 54 + 35) (Ai + A 2 ), 
St (2<?4 + ff 5 ) (Ai + A 2 + A 3 + A 4 ) , 



(111) 

(112) 
(113) 
(114) 
(115) 

(116) 

(117) 
(118) 



(119) 
(120) 



A< B > = 
\C0 - 



-Ai (ff 4 - 2g 5 ) (Ai - A 2 ) , 

-Ai (.94 - 2. 95 ) (Ax - A 2 + A 3 - A 4 ) 



(121) 
(122) 
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\ { 2 ] = 4z (g 4 + 2g 5 ) (A, + A 2 ) , 

A 2 6) = 4z (. 9 4 + 2g 5 ) (A 1 + A 2 + A 3 + A 4 ) , 



(123) 
(124) 



M, 



(123) 



2ao 2c 
2i | &o ao + do c 
2b a 2d 



(125) 



a = [(d + 3)sf 4 + (3d - l).g 5 ] Ai + [(d 2 + 4d - l).g 4 + (3d 2 + 2d - 3)<? 5 ]A 2 

+(d- l)(.g 4 + 3.g 5 )A 4 , 
fe = Vd- l{(.g 4 + 3g 5 )Ai + [(d + 5).g 4 + (3d + 5).g 5 ]A 2 + 5(g 4 + ,g 5 )A 4 } , 
c = Vd - l{(g 4 + 3g 5 )Ai + [(d + 5).g 4 + (3d + 5).g 5 ]A 2 + 2(g 4 + 3.g 5 )A 3 

+ [(d + 7).g 4 + (3d + ll).g 5 ]A 4 + (,g 4 + 3.g 5 )A 5 } , 
do = 5(^4 + ff B )Ai + [(d + 9).g 4 + (3d + 7)g B ] A 2 + 10(.g 4 + .g 5 )A 3 

+ [(d + 19).g 4 + (3d + 17)g B ] A 4 + 5( 34 + g 5 )A 5 , 



(126) 
(127) 

(128) 

(129) 



M. 



(123) 



2a 2 2c 2 
2i | b 2 a 2 + d 2 c 2 
26 2 2d 2 



(130) 



a 2 
b 2 
c 2 



(d 54 + 2 35 ) Ai + [(d 2 + d - l).g 4 + 2dg 5 ] A 2 + (d - l)g 4 A 4 , 
Vd - l{.g 4 Ai + [(d + 2) 34 + 2g 5 ] A 2 + 2(.g 4 + g 5 ) A 4 } , 

Vd - l{ 54 Ai + [(d + 2) 54 + 2.g 5 ] A 2 + 2.g 4 A 3 + [(d + 4) ff4 + 2g 5 ] A 4 + ,g 4 A 5 } . 



d 2 = 2( 34 + g 5 )(Ai + 2A 3 + A 5 ) + [(d + 3).g 4 + 4.g 5 ]A 2 + [(d + 7)g 4 + 8g 5 ]A 4 . 



(131) 
(132) 
(133) 
(134) 



All other b\^} are zero. 



We still have to diagonalize Mq 123 ^ and M 2 123 \ Before explicitly doing so, let us finish the formal development. 
Assume Mq 123 '' to have eigenvectors vqo, v 0+ , vo- with corresponding eigenvalues Aoo, A +, A -. In the same way, 
assume M 2 123 ^ to have eigenvectors w 2 o, v 2+ , v 2 - with corresponding eigenvalues A 20 , A 2+ , A 2 _. The v k and v 2 k arc 
then linear combinations of the u\, u 2 , u 3 . 
If we write the Vj as 



Vq = ^ c ok vok + c\'u 4 + c 6) u 6 . 

fc=0,± 

V\ = c[ 5) u 5 + c^uj , 

V 2 = ^ C 2kV2k + C 2 4) U 4 + 4 6) «6 : 
fe=0,± 



(135) 

(136) 
(137) 



we can write 



B (L) 



^2 c ok^ok v ok + c 4) A 4)L m 4 + c 6) A 6)i u 6 , 



B 



(L) 



fe=0,± 



r (7) A (7)i 7,, 



- B 2 L) = X! C 2fe^2fe 1; 2fe + C. 
fe=0,± 



2 A 2 U4 



c 2 6 ^ \ 2 ^ L uq . 



(138) 

(139) 
(140) 



The final result is then 
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A 20 = 2z(a 2 + d 2 ) , (158) 

A 2 ± = 2i (02 +d 2 ± v / (a2-d2) 2 +4fe 2 c 2 ) . (159) 

Keeping again only quartically divergent terms, they become 

A 20 = 2i[(d 2 + 3d + 4).g 4 + 2(d + 4).g 5 ]AM + 0(A 2 ) , (160) 

A 2+ = 4i(d + 2)[(d + l)g 4 + 2g 5 ] A< 4 ) + 0(A 2 ) , (161) 

A 2 - = 8i( ff4 + 2 55 )A< 4 ) + 0(A 2 ) . (162) 

Now let us consider the rest of the eigenvalues A} . Keeping only quartically divergent contributions, we get 

A< 4) = 16* ( 5 4 + 95) A (4) + 0{A 2 ) , (163) 

A<, 6) = 16* (54 + .95) AW + 0{A 2 ) , (164) 



\ 2 = Si (.g 4 + 2 ff5 )A (4) +0(A 2 ). 
X ( 2 6) =81 ( gi + 2g 5 )A^ +0(A 2 ), 



(165) 
(166) 



In Aj and Aj/ J the quartically divergent contributions cancel. Combining (121) with fl37| ) and (t22h with ( |38| ) and 

(7) 

keeping also quadratically divergent contributions and the leading imaginary part of X\ gives 



A 



(5) 



A 



(7) 



(.94 - 2 55 ) 



(.94 - 2g 5 ) 



A 2 



2{^) 2 m 2 WQ 1 - A 2 V /A 2 W 



•C(A U 



A 2 



8(47r) 2 m 



4 

wo 



4m„ 



HA 2 W /A 2 V ) ^7^(fc 2 +4m 2 W0 )(fc 2 -4m 2 W0 ) 3 / 2 N 
°1-A 2 V /A 2 V + 



3%/PAy 



+0(0,-2), 

where we have used again the shorthand notation (|39| 

divergent are A^ 5 ' and \± ■ In fact, with finit< 
regularization is used merely for convenience 



(167) 



(168) 



Notice that the only eigenvalues that are not quartically 
divergent are A^ 5 ' and A^ 7 '. In fact, with finite Ay and Aw the corresponding integrals are convergent and dimensional 



D. Resonance 



Combining these results with (141)- (144) as well as (146) gives 



S = 



-2i(g 4 - 2g 5 )siu 5 



1 + (54 - 2g s ) 



Al ln(A 2 „/A^) 
2 (4^m^ 1-A V /A^ 



•0(A°) 

-2i(g 4 - 2g 5 )s 1 u 7 



A 2 

1 - (94 - 2# 5 ) 8(4^2^ 



i-2 4 m 2 ln(A^/A 2 v ) CT (fc 2 +4m^ )(fc 2 -4^)3/2 



-0(A" 



0(0,-2) 



(169) 



Under the assumption that the anomalous couplings dominate the gauge coup ling and the addition al as sumption that 
\g 4 — 2g$\ 8(47r) 2 m 4 ,, /(A 2 / m 2 ) [m r is the resonance mass; see eq uati on ( |173| ) below], the terms in (169) are in leading 
order independent of 54 and 55. Modeling the second term in (169) with a Breit-Wigner shape (and consequently 
neglecting the fc 2 -dependence of the width term) gives 
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X r ( SacSbd - S ad Sbc \ ( 9a 1 9% ~ 9aS9% \ 
*, = ^ 2 ' \ 5 I 



+ 



lC 2 



k 2 — m 2 . + im r T r 



S ac 5bd - S ad 6 bc \ ( g^kpks + gfakgkj - gZghpk-t - gf y k a k s 

2k 2 



(170) 



with 



Ci = 4mL ( Hk 2 w /k 2 v ) 



A 2 V \l-A 2 v /A 2 



w 



(171) 



C 2 = 



16(4tt) 



2 mt 



A 2 



(172) 



m r = 2m W0 A 



■0(A" 



and 



7r(m 2 + 4m 2 , )(m 2 - 4m 2 ™) 3 / 2 



Sm^A 2 , 



+ 0(A~ 4 ) . 



(173) 



(174) 



Note that we can trivially drop the superscript "<r" in the C2 term in (170). 

The direct calculation of the resonant part of the bubble sum ( |l7| ) gives the mass ( |173| ) and the width (174) for the 
resonance. As a consistency check for the width, we h ave computed the decay rate of the resonance independently. 
For k 2 « m 2 , the resonant part of the bubble sum (17C) can be written as 



n 



9 



ftU 



k 2 — m 2 + im r T r 



+ Xk^k v \e cdn [(-k) 1 g Sv - {-k)&g lv \ 



with 



A 



2(47r)m 2 



Aym r 

and arbitrary X. That is, the bubble sum can be decomposed into a spin-1 propagator part 

g^v 



and two derivative couplings 



— — : — h Xk^k v 
k z — mi. + im r L r 



^abc(k a gpj — kpg ai 



(175) 



(176) 



(177) 



(178) 



This coupling can be used to compute the decay rate into two W bosons. The result of this standard calculation 
coincides with (174). 
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IV. DISCUSSION 



The first thing one notices about the contributions to the bubble sum is that in most channels quartic divergences 
are present. In these cases we did not calculate the subleading terms. In case quartic divergences are present the 
subleading terms could be of the form A 2 fc 2 . If such terms are absent the whole interaction is suppressed by 1/A 4 
and the channel is of no phenomcnological interest. If such terms are present resonances are present at a scale 0(A 2 ). 
These are out of the reach of present colliders and possibly also of the LHC. Much more interesting is the 1=1, 
,7=1 channel where we find a low- lying resonance in the term 117, which corresponds to longitudinal vector boson 
scattering. Depending on the ratio Ayy/Ay one finds the resonance below or above the vector boson threshold. On 
physical grounds we expect Ay to be the smaller, as it is directly related to the Goldstone boson sector of the theory, 
where the strong interactions are supposed to take place. In that case the resonance always lies above the two-W- 
threshold. Because the interactions of the transversal vector bosons are suppressed by the gauge coupling, a reasonable 
assumption would be Ay w gA\y- This corresponds to m r w 200GeV. A recent comparison with the LEP-100 data 
gives a limit Ay > 490GeV 0. For m r = 200GeV, this gives T r < 12GeV. The fact that the resonance can be at 
such low energy is somewhat surprising, given our experience with chiral perturbation theory in pion physics. 

To study the connection with pion physics we make the substitutions 174 = g 4 £4, .95 = g 4 e$ and m wo = gf^/2. In 
the resulting Lagrangian one takes g — > 0, with £4, £5 and /„• fixed. This way one finds the standard nonlinear sigma 
model with two higher derivative terms. We define g p = 4(f4 — 2£ 5 ). For didactical purposes we keep in the chiral 
perturbation theory here the tree level terms, the imaginary part of the ordinary chiral loop and the contribution of 
the first loop we calculated. One finds for the I = 1, J = 1 amplitude, 



After unitarizing the amplitude with the [1, 1] Pade approximant one finds the p resonance with a width given by 



We had to ignore the s term in the imaginary part of an because it is of too high an order in chiral perturbation 
theory and other terms of this order have been ignored. The presence of the extra loops that we calculated results 
therefore in a correction to the KSRF || relation T p = m^/ (96-Trf 2 ). Of course, this is not the full story for chiral 
perturbation theory, but it shows that the effects we calculated are only a small correction to pion physics. The reason 
we get a large effect in our calculation is that we assume that the interactions are dominated by the g p interaction. 
In the chiral limit this is not possible, because at low enough energy the lowest order term in s always dominates. 
Since within the standard model there is the W threshold to consider, it makes sense to say that the anomalous term 
dominates. It means \g p \s 2 /v 2 3> s for s > 4to 2 , . Within the standard model this translates into A\g4 — 2g 5 \ 3> g 2 . 
This condition therefore quantifies what we mean by strong anomalous couplings. 

The model we discussed fits in naturally in the class of models that give rise to resonances in W physics due 
to strong interactions. This class of models is collectively known as the BESS model, from breaking electroweak 
symmetry strongly. A recent review of this class of models is jnj. Nonetheless, as the discussion above points out, the 
model is subtly different from the models in the literature, due to the dominance of the anomalous couplings. This 
would make the model into a mere curiosity were it not for a very simple class of models giving rise to precisely the 
required structure. This is the class of the strongly interacting singlet Higgs (SISH) models M. These models, which 
contain beyond the standard model only extra scalar singlet particles, are actually the simplest possible renormalizable 
extensions of the standard model. Because the singlets couple only to the Higgs boson, they do not change the 
phenomenology at LEP-100 at the one-loop level. Two-loop effects are too small to be significant. In order to perform 
a precise phenomenology of the model, SUr(2) breaking effects should be taken into account. The model satisfies 
the LEP-100 limits on extra Z bosons, because both the coupling to leptons and the mixing with the Z boson are 
suppressed by an electroweak loop. The statistics at the Tevatron is too small to see the resonance. For the planned 
high-energy colliders the phenomenology should be straightforward, the resonances being produced via vector boson 
fusion. For LEP-200 the situation is somewhat subtle, as the resonance does not couple to the incoming electrons 
directly. Strong form factor effects could play a role. The precise phenomenology will be left to future work. 

Given the fact that the couplings are strong, one can question how generic the results are. In principle, higher-order 
terms in the chiral perturbation theory could be important. For the formation of the resonance via a bubble sum, 




(179) 




(180) 
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only four-point vertices contribute. The presence of higher order terms will therefore not effect the structure of the 
calculation very much. Basically, we expect a form factor for 54 — 2g 5 which, when inserted into the graphs, makes 



the explicit dependence on Ay more complicated. However, the term (g 4 — 2g 5 )Ay in formula (169) is essentially 



v 

determined by power counting; so one would expect corrections of the form A 2 — > A 2 + C(4m 2 v ). Also, the precise 
formula for the resonance mass as a function of Aw, Ay could become more complicated. However that will not 
change the qualitative picture of m r ~ 200GeV, with a coupling of order m r /Ay. 

Finally one can ask if the model is consistent with the LEP precision data. The limits derived on 54, g§ in HQ follow 
from simple perturbation theory and depend crucially on the behavior in the hypercharge sector of the theory, as they 
come from the limits on the p parameter. For the strongly interacting case as described here they are not expected 
to be a good estimate and could possibly even be an order of magnitude wrong. The model is similar in appearance 
to the standard model in the large Higgs boson mass limit, which is disfavored by the LEP1 data. When one simply 
removes the Higgs boson from the standard model, the model becomes nonrenormalizable, but the radiative effects 
grow only logarithmically with the cutoff. The question is whether this scenario is ruled out by the LEP1 precision 
data. The LEP1 data appear to be in agreement with the standard model, with a preferred low Higgs boson mass. 
One is sensitive to the Higgs boson mass in three parameters, known as S, T, U or ei, £2, £3. They receive corrections 
of the form g 2 [hi(mH /inw) + const], where the constants are of order 1. The logarithmic enhancement is universal 
and would also appear in models without a Higgs boson as ln(A), where A is the cutoff where new interactions should 
appear. Only when one can determine the three different constants can one say one has established the standard 
model. At present, the data do not suffice to do this to sufficiently high precision. In practice, one can compensate a 
change in the Higgs boson mass in the formulas with extra contributions to the S, T and U parameters. As such terms 
are generated by the contributions of formula (||), there is enough freedom to fit the data; see |?j for a discussion. 
Whether a model with a low-lying resonance would actually improve the fit to the data depends on the couplings to 
the fcrmions and to the hypercharge. 
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APPENDIX A: ONE-LOOP RESULTS 

Here we list results for the necessary basic one-loop integrals and for the composite one-loop quantities Ai through 
A 5 . 

Let d be the spacetime dimension, d — 4 — 2e, let \i be given by ln47r/i 2 — 7 B = ln/x 2 with the Euler-Mascheroni 
constant 7^ and use the abbreviation J given in ( p3| ) . 



The only integrals we need are 



1 

p p 2 — to 2 (47r) 2 \e /2 2 



n™ 2 )= I -^-, = ^(- + l-ln^)+0(e) (Al) 



and 



J(fc 2 ;m 2 ,m 2 )= / - (A2) 

J P [(P + k) 2 — rn? a + ie] (p 2 - m 2 + is) 



Define 

D = fc 4 + m 4 + to 4 - 2fc 2 m 2 - 2fc 2 m 2 - 2m 2 m 2 . (A3) 
For k 2 > we get the following results in a straightforward calculation. 
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1. I (k 2 ; mi, m 2 b ) for D < 

D < is equivalent to the statement that none of fc 2 , m 2 ., to 2 is larger than the sum of the other two. 



I(k \ m a ,m b ) 



(4tt) 2 



k 2 + to. 



2 2 2 

- — mi . to; 



2k 2 
arctan 



In- 



k 2 + m? — to. 



2fc 2 



fc 2 + to 2 — to 2 



arctan - 



k 2 + m? — 



0(e) 



I(k 2 ;m 2 a ,m 2 b ) 



(4tt) 2 
+0(e). 



2. /(fc 2 ;m 2 ,m 2 ) for D > with fc 2 < |m 2 - m 2 

i 2 - fc 2 ln ln '"-q T 



to 2 — m 2 — fc 2 , to 2 to 2 — m 2 — k 2 to? \AD . m 2 + m? — fc 2 + \AD 



2fc 2 



■ln-i 



2fc 2 



& 2 TO 2 + TO 2 - fc 2 - \/£> 



Equation (A5) can also be written as 



Til 2 2 2\ 



(4tt) 5 



or 



/(fc ;m a ,m b 



(4tt) 5 



1 2 w 2 - m 2 a - fc 2 + VP ^ m| w 2 - to 2 - fc 2 + yg ^ to| 

~ 2 2fc 2 u 2 



fc 2 



in 



2fc 2 




TO 2 + TO 2 — fc 2 + 




1 2/I 2 




2 2 7 2 

TO b — m a — fc — 


^1 


2fc 2 




TO 2 + TO 2 — fc 2 — 





■0(e) 



ln^ 



2fc 2 



ln^ 

/i 2 



fc 2 



2/Z 2 



0(e) 



3. /(fc 2 ;m 2 ,m 2 ) for D > with fc 2 > |m 2 - m 2 j 



Here 



Til 2 2 2 

/(fc ;m a ,m b 
1 




167rfc 2 
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4. I(fc 2 ;m 2 ,ra 2 ) 



If D = k 2 (k 2 - Am 2 ) < 0, i.e., for k 2 < Am 2 , we get from (|A4| ) 

1 



I(k z ;m 2 ,m z ) 



1 1 2 1 5 U + 2 - ln^ - 2*/^ - 1 arctan 



(4tt) 2 



1 



•O(e) 



For fc 2 <C to 2 , we can expand in powers of k 2 /to 2 to get 
I(k 2 ; to 2 , to 2 ) 



(47T) 



1 , to 2 1 / /fc 2 

~ ~ ln ~ + ft ~ 2 

e /i z d y to 



60 V 



■O e 



If D = k 2 (k 2 - Am 2 ) > 0, i.e., for k 2 > Am 2 , ( |A8[ ) becomes 



I(k 2 ;m 2 ,m 2 ) 



(Any 



1 , to 

- + 2 - ln — 

e /i z 



1 - 



4m 2 



l - 



In- 



-a/1- 



4 m- 



2m 2 



1 - 



4 m 2 



4m 2 



16tt 



(A9) 



(A10) 



•0(e). (AH) 



5. I(fc 2 ;m 2 ,0) for fc 2 < m 2 



Now -D = |fc 2 — to 2 | > 0. Let us assume k 2 < to 2 . Then we get from (A6) or (A7) 



I(k 2 ;0,m 2 ) = l(k 2 ;m 2 ,0) 



(An) 2 
i 

jAnj 2 



- + 2- 
e 



k 2 



ki- 



rn TO 



1 to 2 / k 2 
- + 2 + —, r In 1 =■ 



/fc 2 



-In- 



0(e). 



For fc 2 <C to 2 , we can expand in powers of k 2 /m 2 to get 
7(/c 2 ;0,to 2 ) = 7(fc 2 ;TO 2 ,0) 



(An) 



1 , TO 2 

- + 1 - In — 



1 ftf_ 

2 Ito 2 



1 fk 2s 1 



6 V to 2 



1 /fc 2x ; 



1 /fc 2 ^ 4 ' 



12 \ to 2 7 20 V to 2 



(A12) 



(A13) 



6. /(fc 2 ; m 2 , m. 2 ) for fc 2 , m 2 -C m 2 



If fc , to 2 <§; to 2 , we can expand (A6) or (A7) in negative powers of m 2 , to get 



I(k 2 ;m 2 a ,m 2 b ) 



(An) 2 



_ + i _ i n _£ 

e /i z 



ifc 2 + to 2 In =f fc 2 (±fc 2 + |to 2 ) + to 2 (A: 2 + to 2 ) ln =8 



+ - 



k 2 (±k 4 + \k 2 m 2 a + §to 4 ) + TO 2 (fc 4 + 3fc 2 TO 2 + to 4 ) ln =j 



+ - 



fc 2 (^fc 6 + f§fc 4 m 2 + ffc 2 ™ 4 + |m«) + TO 2 (fc 6 + 6fc 4 rn 2 + 6fc 2 TO 4 + m 6 a ) ln ■ 



+0(e,TO h - 10 ln to 2 ) 



(A14) 
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7. I(k 2 ; ml, ml) for k 2 -C m^mj 



Here 



I {k 2 ; m 2 a , mf) 



(47T) 



e 



+G(fc 8 ,e). 



TOf 



2(m 2 



m 2 ) 2 



("»« 



m\ + lOmlml + m^ m l m li m l + m l) ^ m l 



6(to 2 



2 - to 2 ) 4 



{ml 



2 ~ ml) 5 



29to 4 to? + 29m 2 ,mf 



12(m 2 



if 



4 _ m a m K m a + im 2 a m 2 b + ml) ln ™a ) k 6 



(TO 2 - TO 2 ) 



(A15) 



8. Calculation of Ai, A 2 , A 3 , A 4 , A 5 

To compute the quantities Ai through A5 defined in (p4|), we first need some preliminaries. Let / in the following 
be some function of p 2 and (p + k) 2 . 
Contracting 

jpaf = A ^ (A16) 

with k a gives 

A = /"(p • k)f . (A17) 

Contracting 

Jp a p f = A' 9a0 + B'^ (A18) 
with <7 Qj a and k a kp and solving the resulting equations for A' and B' we get 

^^/(V-^)/, (A19) 

B'=«hJ p (-P 2 + d^)f. (A20) 

Contracting 

/ PaP0P 1 f = A + B — — — (A21) 

Jp K K 

with fc Q <7^ 7 and fc a fc^fc 7 and solving the resulting equations for A" and B" , we get 

A " = ^lj (P 2 (p-k)-^)f, (A22) 
B" = -± J J^-3p 2 (p-k) + (d + 2)^)f. (A23) 

20 



Contracting 



gapk-fks+gctykgks+gaSkpk^ 

r All// i i \ i r,m +g-,sk a kp+g l3s k c ,k- l +gp- l k a k s _ ;// k a knk~,k s 

PaPpP-fPsf = A"'(g a pg^ 5 + g ai g 0S + g a &g 01 ) + B p h C" — 



(A24) 



with g a /3g 7 s, k a k g 7 s and k a kpkyks and solving the resulting equations for A'" , B"' and C" we get 

A'" = VrmW)j p (p 4 + ^)f> (A25) 

B'" = TS^ml (-P 4 + (d + - (d + 2)^) / , (A26) 

C "' = (d -i)Wi) / ( 3p4 " 6 ^ + 2 ) £ ^ + ( d + 2 )( d + 4 )t) / • (A27) 

Using the above results, some trivial but lengthy algebra on the integrand, the properties of dimensional regular- 
ization, the results for one-loop integrals in the preceding parts of the appendix, as well as 

(P • k) 2n f(p 2 ) = c n (k 2 ) n [ (p 2 ) n f(p 2 ) (A28) 



with 



i.e. 



2n — 1 

co = 1 , c n = — — -c n _i , (A29) 

a + 2n — 2 



_ r(n + l/2)r(d/2) 
" " r(n + d/2)r(l/2) ' 1 j 

we can compute Ai through A5. The results are given in the main text in eqs. (pq) through (Ba). 
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